Nonlinear Schrodinger equation with chaotic, random, and nonperiodic nonlinearity 
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In this paper we deal with a nonlinear Schrodinger equation with chaotic, random, and nonperiodic cubic 
nonlinearity. Our goal is to study the soliton evolution, with the strength of the nonlinearity perturbed in the 
space and time coordinates and to check its robustness under these conditions. Comparing with a real system, 
the perturbation can be related to, e.g., impurities in crystalline structures, or coupling to a thermal reservoir 
which, on the average, enhances the nonlinearity. We also discuss the relevance of such random perturbations 
to the dynamics of Bose-Eisntein Condensates and their collective excitations and transport. 
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The nonlinear Schrodinger equation (NLSE) is the math- 
ematical vehicle that describes the evolution of solitonic so- 
lutions for different nonlinear systems, such as, fiber optics 
HI], bulk medium and photonic crystals 0], Langmuir waves 
in plasmas H], wave function of Bose-Einstein condensates 
(BECs) Jl, and others. 

A special case involving the NLSE consists in variable co- 
efficients modulated in the spatial and/or temporal coordi- 
nates. The control of these coefficients allow us to obtain new 
distinct solutions. In this context, JH @] have recently pro- 
posed a treatment of BECs using similarity transformations to 
construct explicit nontrivial solutions of the cubic and cubic- 
quintic NLSE with potentials and nonlinearities depending 
both on time and on the spatial coordinates. Also, thermal 
effects on nonlinearities can change its form, presenting new 
solutions 1 7]. In BECs, the s-wave scattering length of inter- 
atomic collisions determines the strength of the nonlinearity 
coefficient |4|], and it can be controlled using the Feshbach 
resonance (FR) [8] via external magnetic Si or optical 

in] 

fields. The FR mechanism allows for a practical means 
to manipulate the nonlinearity 111 211 . 

Although the control of the nonlinearity has been very ef- 
fective, noises can appear in the system management or added 
to it. In this way, the perturbations can change the nonlinear- 
ity, and thus influencing possible changes in the noise-free so- 
lutions. The inclusion of spatial random potential in the study 
of BEC dynamics has been proposed in Refs. II131 fl4l fl5l fl6n 
and recently was tested experimentally ifTvh . When random- 
ness is introduced in the BEC dynamics, through an optical 
speckle, one may be able to study Anderson localization in 
the context of BEC and superfluidity H 1 311 . It was clearly 
demonstrated in these papers that in the presence of disorder 
the condensate's expansion in ID waveguides is inhibited, and 
the collective dipole and quadrupole oscillations are strongly 
damped. In a way, this is similar to the damping of collective 
states in nuclei and metal clusters [18], where the randomness 
is internal. 

A natural question arises as to what would be the conse- 
quence of having the disorder present directly in the nonlin- 
earity term? For BEC, this implies a point two-body inter- 



action (t-matrix) with a random component. Would this add 
or remove some of the effects of the speckle potential? In 
this connection, [19] have recently proposed an NLSE in the 
presence of random nonlinearity. Specifically, these authors 
consider the effects of random time modulation of the non- 
linearity coefficient on the dynamics of solitary waves in the 
NLSE. On the other hand, to our knowledge, chaotic perturba- 
tions in the nonlinearity term in the NLSE have not been fully 
considered yet, though it is a common knowledge that several 
physical systems do exhibit chaotic behavior. 

Classically chaotic systems appear to behave as random 
systems. Tiny differences in the initial state of the system can 
lead to enormous differences in the final state even over fairly 
small time scales W20I1 . This happens even though these sys- 
tems are deterministic, meaning that their future dynamics are 
fully determined by their initial conditions with no random 
elements involved. This behavior is known as deterministic 
chaos. 

The dynamical systems theory (DST) is an area whose in- 
terest lies mainly in nonlinear phenomena, the source of clas- 
sical chaos. DST groups use several concepts to the study of 
chaos, such as Lyapunov exponents, fractal dimension, bifur- 
cation, and symbolic dynamics among other elements lEoll . 
Recently, other approaches have been considered, such as in- 
formation dynamics and entropic chaos degree 112 ill . For ex- 
ample, given Co, C n can be the n-th iterate of the quadratic 
functions: C n (/x) = C^ L _ 1 + \i\ sine functions: C n (/i) = 
/isin(C n — 1); logistic functions C n (fi) = /iC n _i(l — C n _i); 
exponential functions: C n (fi) = /iexp(C n _i); doubling 
function defined on the interval [0, 1): C n = 2C n _i mod 1, 
and so on, fi being a parameter. It is worth recalling that all the 
functions in the above list are familiar to researchers in DST. 
For example, for some values of /i, it is known that some of 
these functions can behave in quite a chaotic manner [20] . In 
what follows we make a distinction between chaos and ran- 
domness, though both concepts indicate disorder. 

The major thrust of our paper is to verify the influence of the 
different types of perturbations in the nonlinearity of a system 
governed by NLSE. In this sense, we investigate the chaotic, 
random, and nonperiodic nonlinearity perturbation. We know 



2 



that these perturbations are different, however, to what extent 
the overall effect is universal (independent on the details of 
the random perturbations), and how can it modify a solitonic 
solution? We purports to supply some answers to the above. 

Differently from Ref. W19I1 . which considers a random time 
modulation on a certain point (generating a Gaussian distri- 
bution), here we consider a constant background nonlinearity 
perturbed by a random function that interferes in both spatial 
and temporal coordinates. Surprisingly, some of the solutions 
found here can move depending of the amount of perturbed 
points in the nonlinearity. This fact is similar to those studied 
by the thermal effects on the nonlinearity |7|]. 

Firstly we consider the NLSE given by 

iipt = -^ xx + g(x, t) \^\ 2 ip, (1) 

where ip = i/j(x, t), x, and t dimensionless, and g(x,t) is the 
function that describe the nonlinearity of the system. Here we 
consider 

g(x,t) = G(l + a(x,t)), (2) 

where G and Ga(x, t) are the nonlinear parameter and the co- 
efficient generated by a chaotic, random, or nonperiodic gen- 
erator, respectively. Eq. (OQ) describes, e.g., the density of par- 
ticles in a Bose-Einstein condensates when it is free of exter- 
nal potentials, in a configuration type cigar- shaped; the spatial 
pulse propagation in bulk crystals that present Kerr-effect in 
the nonlinearity or temporal pulse propagation in nonlinear 
optical fiber; etc. 

We investigate the evolution of the solution for the NLSE 
via numerical simulations, based on the split- step finite dif- 
ferences (SSFD) method with a time-step and space-step size 
of At = 0.0001 and Ax = 0.01, respectively; we use N to 
represent the number of points in space. The core of SSFD is 
based on the Crank-Nicolson algorithm 12211 . To control these 
numerical simulations we looked for the conserved quantity 
given by (power) 

N 

P = Y J \^t)\ 2 . (3) 

x=l 

To calculate the error we use the comparative form 

1 N 

Er = ^^(mxM 2 -m^tf)\ 2 ), (4) 

x=l 

where tf is the final time of the evolution. The equation above 
represents a mean distance between the input and output state. 
When we take G = — 2 and a(x,t) = in ©, we can write 

ip = e ifit sech(x) (5) 

as solution of (Q]), with [1=1. This solution will be taken as 
initial condition for our simulations of the Eq. (Q]). 

For our numerical simulations we consider the chaotic per- 
turbation in the nonlinearity given by the logistic function, the 




FIG. 1: Plots of the perturbed nonlinearity for the (a) chaotic, (b) 
random, (c) and nonperiodic functions. 



random perturbation is generated by random algorithm simu- 
lator, and the nonperiodic perturbations are generated by the 
function a(cos(5x)/2 ± cos(Vbx)/2), where a assumes the 
values of the perturbation. In Figs. Hkttb, andQJ we display a 
form of the chaotic, random, nonperiodic nonlinearity as func- 
tions of space for a generic time, respectively. Experimentally, 
this perturbation in the nonlinearity can be constructed, e.g., 
in a crystal with impurities which are altered chaotically, ran- 
domly, or non-periodically. It remains to be seen what the 
presence of Ga(x, t) implies, though one would guess that it 
amounts to taking into account, within the mean field, Gross- 
Pitaevskii, description, the effects of the many-body correla- 
tions. 

We use © as input state in (Q]) to verify its evolution in the 
presence of the above mentioned perturbations. In Fig. [2^ we 
consider the chaotic perturbation via a(x,t) between ±0.001 
and ±0.1 (10%). The chaotic perturbation is obtained consid- 
ering 4000 points affected within the interval — 20 < x < 20 
in space, that changes the spatial profile of the nonlinearity. 
These perturbed points are changed by a new function after a 
time t = 20, and so we will have 100 temporal points affected 
in the interval < t < 200. Fig.[5J) displays the soliton am- 
plitude (height) of the solution (|^| 2 ) at the position x = 0. 
Note that it becomes vanishingly small after t = 60 owing to 
its motion. However, the soliton is stable in this range of per- 
turbation. We calculated the error in the power of2.41xl0 -4 
and the comparative error of E r = 9.39. 

On the other hand, when we consider a chaotic perturba- 
tion of 50% in the value of the amplitude, we found that the 
soliton practically disappears. This case is shown in Fig. [3] 
In the Figs. [3^ and[3j) we display the |?/;| 2 and the height in 
the position x = 0, respectively. The error in the power was 
2.38 x 10 -4 with a mean distance E r = 9.76. In this case the 
mean distance is due to the moving and the vanishing pattern, 
i.e., the output state can not be at the same position as that of 
the input state. 

Now, when we consider a random perturbation of the non- 
linearity, obtained here by an algorithm of random number 
generation, we note a different behavior of the solutions. Here 
the soliton remains stable even for 50% of perturbation, differ- 
ently from the case of the chaotic perturbation. In Fig.Hk we 
plot the soliton solution considering 10% of random pertur- 
bation in the nonlinearity using the same arguments presented 
for the chaotic case, i.e., 4000 affected points in space versus 
100 temporal points into the range shown. Fig. 0J) displays 
the amplitude of the soliton at x = 0. This perturbation is 
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FIG. 2: Plots of the soliton evolution \tp(x, t)\ 2 with 10% of chaotic 
perturbation in the nonlinearity. In (a) is displayed the solution and 
its profile (top panel) and (b) its height at position x = 0. 
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FIG. 3: Plots of the soliton evolution \ip(x, t)\ 2 with 50% of chaotic 
perturbation in the nonlinearity (see text for details). In (a) is dis- 
played the solution and its profile (top panel) and (b) its height at 
position x = 0. 



responsible for moving the soliton. The error in the power is 
of 4.88 x 10 -6 with the mean distance between the input and 
output states of 10 _1 . For 50% of random perturbation the os- 
cillation of the soliton is more evidenced when compared with 
the case of 10% of random perturbation. Figs.Ek and0) show 
the \i(j\ 2 and the amplitude for x = 0, respectively. The ro- 
bustness is guaranteed, differently form the chaotic case. The 
errors in power and the mean distance are 1.17 x 10 ~ 5 and 
10 _1 , respectively. 

To conclude our study we investigate the effects of the non- 
periodic perturbation in the nonlinearity. From Fig. \Tj^ one 
observes that this perturbation seems the most well-behaved 
compared to the other two types. This fact is reflected in the 
behavior of the solution which remains practically with the 
same form as that of the input state. This occurs even when 
it suffers 50% of nonperiodic perturbation, as can be seen in 
Fig- El The error in the power is of 2.83 x 10 -4 and the mean 
distance is of 2.88. 

With the results presented here, we verify that, under simi- 
lar conditions, chaotic, random, and nonperiodic perturbations 
in the nonlinearity can present distinct features, and some- 
times results in vanishing solitons, as verified when the system 
suffers chaotic perturbations. 

In summary, in the present work we have studied the ef- 
fects of chaotic, random, and nonperiodic perturbations in the 
nonlinearity on the soliton evolution via NLSE. We consid- 
ered cubic nonlinearity with strength perturbed chaotically, 
randomly, or nonperiodically. In the chaotic case we found 
that moving solitons can be destroyed when they are perturbed 
with 50% in the intensity of the nonlinearity. On the other 
hand, when the system engenders random perturbation this 
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FIG. 4: Plots of the soliton evolution \ip(x,t)\ 2 with 10% of ran- 
dom perturbation in the nonlinearity (see text for details). In (a) is 
displayed the solution and its profile (top panel) and (b) its height at 
position x = 0. 
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FIG. 5: Plots of the soliton evolution |^(x,£)| 2 with 50% of ran- 
dom perturbation in the nonlinearity (see text for details). In (a) is 
displayed the solution and its profile (top panel) and (b) its height at 
position x = 0. 



does not occur. The soliton solution remains stable, how- 
ever now it can move. Finally, when we look for the non- 
periodic perturbation we found that it displays robust solu- 
tions with no apparent influence on the solitons. In a way, 
we can say that disorder in the non-linearity may, or may not 
lead to Anderson-type localization, depending on the nature 
of the perturbation. Our results have direct impact on work in 
optical lattices with impurities in the crystal, laser-generated 
randomness in the non-linearity, and many-body effects in the 
dynamics of Bose-Einstein condensates and their collective 
excitations and transport. 
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FIG. 6: Plots of the soliton evolution \i/j(x, t)\ 2 with 50% of nonpe- 
riodic perturbation in the nonlinearity (see text for details). In (a) is 
displayed the solution and its profile (top panel) and (b) its height at 
position x = 0. 
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